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Abstract 

We report results of a search for CP violation in the singly Cabibbo-suppressed 
decays D+ -> R-K+-K+, ^vr+, K*(892)°K+, and 7T 7T + 7r + based on data from the 
charm hadroproduction experiment E791 at Fermilab. We search for a difference in 
the D + and D~ decay rates for each of the final states. No evidence for a difference 
is seen. The decay rate asymmetry parameters (Aqp), defined as the difference in 
the D + and D~ decay rates divided by the sum of the decay rates, are measured to 
be: A C p(KKit) = -0.014 ± 0.029, A C p(<H = -0.028 ± 0.036, A C p(K*(892)K) = 
-0.010 ± 0.050, and A C p(innc) = -0.017 ± 0.042. 

PACS: 11.30.Er; 12.15.Ji; 13.25.Ft; 14.40.Lb 



CP violation can be accommodated in the Standard Model (SM) by a complex 
phase in the Cabibbo-Kobayashi-Maskawa matrix describing the transitions 
between quarks induced by the charged weak interaction. To date, CP viola- 
tion has been observed only in the neutral kaon system. It is being actively 
searched for in B decays, where the effects of CP violation in the SM are 
expected to be large |TJ, and in hyperon decays 0. In contrast to the strange 
and bottom sectors, the SM predictions of CP violation for charm decays are 
much smaller making the charm sector a good place to test the SM and 
to look for evidence of physics beyond the SM [|J. 

CP violation occurs if the decay rate for a particle differs from its CP con- 
jugate decay rate: Y{D — ► /) ^ T(D — > /). Such an asymmetry requires 
the interference of at least two independent amplitudes with non-zero relative 
phase. In D decays, as in K and B decays, there are two classes of CP viola- 
tion: indirect and direct. In the case of indirect CP violation, the asymmetry 
is associated with mixing, which can occur only in neutral meson decays. In 
the case of direct CP violation, final state interactions must induce a strong 
phase shift and both the strong and weak phases of the two amplitudes must 
differ. This can occur in both charged and neutral meson decays. In the SM, 
direct CP-violation asymmetries in D decays are predicted to be largest in 
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singly Cabibbo-suppressed (SCS) decays (at most of the order of 10~ 3 ) and 
non-existent in Cabibbo-favored (CF) and doubly Cabibbo-suppressed (DCS) 
decays ||. 

Current experimental sensitivity to decay rate asymmetries due to direct CP 
violation, defined below, is of the order of 10 _1 . Experimental results on 
searches for CP violation in D° decays come from E691 ||, E687 ||, and 
CLEO |7p. However, limits on CP violation in D + decays come only from 
E687 ||. In this letter, we report on a higher-statistics search for evidence 
of direct CP violation in the SCS decays D + — > K~K + tc + (inclusive), 4>ti + , 
K*(892)°K + , and, for the first time, tt-it + tt + . 

The signature for CP violation in charged D decays is an asymmetry in the 
decay rates: 

_ r(p+->/+)-r(fl-->/-) 

CP T(D+ -> f+) + T(D~ -> /-)' {) 



To the extent that 



'(!>■ - tics) <D~ - fs CS ) 



e(D+ -> K-1T+TT+) e(D~ -> K+ 



7T 7T 



(2) 



where e is our detector efficiency and fscs ls K^K ± ti ± ) (jnr^, K*(892)K ± 1 
or 7r =F 7T ± 7r ± , we can normalize the production rates relative to the CF mode 
D ± —>■ fr =F 7r ± 7r ± . Then, Eq. ([!]) becomes 

V (D+) - V (D-) 

AcP - r,(D+)+r,(D-y (3) 
where 

and iV is the number of observed D decays. Thus, differences in the D ± pro- 
duction rates ||, relative efficiencies, and most other sources of systematic 
errors cancel. 



E791 is a high statistics charm experiment which ran at Fermilab during the 
1991-1992 fixed-target run ||. The experiment combined an extremely fast 
data acquisition system with an open trigger to record the world's largest 
sample of open charm. Over 20 billion events were collected using the Tagged 
Photon Spectrometer with a 500 GeV n~ beam. There were five target foils: 
one 0.5 mm thick platinum foil followed by four 1.6 mm thick diamond foils 
with 15 mm center-to-center separations. The spectrometer included 23 planes 
of silicon microstrip detectors (6 upstream and 17 downstream of the target), 
2 dipole magnets, 10 planes of proportional wire chambers (8 upstream and 
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2 downstream of the target), 35 drift chamber planes, 2 multicell threshold 
Cerenkov counters that provided tt/K separation in the 6-60 GeV/c momen- 
10 , electromagnetic and hadronic calorimeters, and a muon de- 



tum range 
tector. 



All selection criteria (cuts) used to choose candidate SCS decays — > 
K T K ± n ± , 4>n ± , K*(892)K ± , and 7r =F 7r ± 7r ± , with the exception of Cerenkov 
identification requirements, were optimized to maximize S/ \/S + B, where S 
is the number of CF D ± — > K T n ± 7i ± signal events, scaled to the level expected 
for each SCS decay mode, and B is the the number of background events for 
each of the signals studied. To reconstruct the D ± candidates, three-prong 
decay vertices with charge of ±1 were selected. All decay tracks were required 
to travel through at least one magnet and be of good quality. Decay vertices 
were required to be located outside the target foils, and the significance of 
separation from the primary vertex in the beam direction (Az/a^z) was re- 
quired to be at least 11 for KKn, 8 for </>tt, 9 for K*K, and 14 for tttttt. The 
component of the D momentum perpendicular to the line joining the primary 
and secondary vertices was required to be less than 0.35 GeV/c for KKn, (fin, 
and K*K, and 0.20 GeV/c for tttttt. The impact parameter of the D momen- 
tum with respect to the primary vertex was required to be less than 55 fim for 
KKn and K*K, 80 /zm for 0tt, and 40 fim for tttttt. Decay tracks were required 
to pass closer to the secondary than to the primary vertex and the sum of the 
square of their momenta, perpendicular to the D direction, was required to 
be larger than 0.30 (GeV/c) 2 for KKn, K*K, and tttttt, and 0.15 (GeV/c) 2 

for 0TT. 



In the decay D — > K T n n , the kaon was identified on the basis of charge 
alone and no Cerenkov identification cuts were applied. Also, in the decay 
D ± — > (fi>n ± , no Cerenkov identification cuts were applied due to the (fi mass 
selection criteria. In the decays — > K^K^n^ and i^*(892)i^ ± , Cerenkov 
identification cuts improved the signal significance. In these decay modes, 
tracks considered as kaons were required to have kaon signatures in the Cerenkov 
counters. No Cerenkov identification cuts were applied for the decay — > 
tt =f tt ± tt ± . To minimize the systematic uncertainties when extracting Acp, we 
applied the same cuts, with the exception of Cerenkov identification, to the 
SCS decays and to the normalizing CF decay for each of the decay modes 
studied. 



The mean momentum of surviving D^ 1 mesons was 70 GeV/c. Longitudinal 
and transverse position resolutions for the primary vertex were 350 /im and 
6 /im, respectively. For secondary vertices from D decays, the transverse 
resolution was about 9 /im, nearly independent of the D momentum, and 
the longitudinal resolution was about 360 /im at 70 GeV/c and increased by 
30 /im every 10 GeV/c. 
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To determine the yields for each decay mode, we performed simultaneous 
binned maximum likelihood fits of the mass distributions (see Figs. for 
the positive and negative candidates. The spectra were fitted with Gaussian 
signals and linear backgrounds. The widths of the Gaussian functions were 
constrained to be the same. Reflections from — > K^tt tt appeared near 
the Df — > K t K ± tc ± and K*(892)K ± signals, and as a broad shoulder below 
the — ► 7r =F 7r ± 7T ± signal. These reflection regions were excluded from the 
fits. 

The mass spectra for all K^K^'K combinations that survived the optimized 
cuts are shown in Fig. [|. The simultaneous fits of the K~K + tt + and K + K~n~ 
mass spectra yielded 1031 ± 44 and 1265 ± 48 signal events for the D + and 
D~ , respectively. 

To reconstruct the in the decay modes 07r ± with <ft —>■ K + K~ , we required 
candidates to have M(K + K~) within ±6 MeV/c 2 of the <fi mass. Because of 
the real K + K~ background under the <fi signal, we performed <j) mass sideband 
subtraction to determine the true number of candidates decaying to 07r ± . The 
sideband regions, which were chosen to have summed area equal to the 
background in the signal region, were 0.990 < M(K + K~) < 1.000 GeV/c 2 
and 1.040 < M(K + K~) < 1.044 GeV/c 2 . The simultaneous fits of the 0vr+ 
and 4>7r~ mass spectra, shown in Fig. yielded 474 ± 25 and 598 ± 28 signal 
events for the D + and D~, respectively. 

To reconstruct D ± -> K*(892)iT ± with iT(892)° -> R-n + and K*(892)° -> 
K + 7i~, we required candidates to have MlK^ir^) within ±45 MeV/c 2 of the 
if* (892) mass. Here, we also performed if* (892) mass sideband subtraction 
to determine the true number of candidates decaying to K*(892)K ± . The 
if* (892) sideband regions were 0.680 < M( J fsT T 7r ± ) < 0.750 GeV/c 2 and 
1.040 < M( J fsT T 7T ± ) < 1.070 GeV/c 2 . The simultaneous fits of the K*(892)°K+ 
and K*(892)°K- mass spectra, shown in Fig. |, yielded 239 ± 18 and 291 ± 19 
signal events for the D + and D~ , respectively. 

In the study of the decay — > 7r T 7r ± 7r ± , we investigated possible reflec- 
tions from D^ 1 — > K t 7t ± ti ± and K T K ± Tt ± , and from random combinations of 
D°(D ) -> ^vr ± , tt^tt* and K^K± with a 7r . The only significant reflection 
was from — >• i^ =F 7r ± 7r ± . We excluded this region from the fits, as shown 
in Fig. The 7r _ 7r + 7r + and 7r + 7r~7r - mass spectra were simultaneously fit us- 
ing four Gaussian functions for the D ± and Df 7r T 7r ± 7r ± signals. The fits 
yielded 697 ± 42 and 851 ± 48 signal events for the D + and D~, respectively. 

We determined the yields of the CF normalizing decay mode — > K^^it^ 1 
that survived the optimized cuts for each of the SCS decays D ± — > K T K ± n ± , 
(p^, K*(892)K ± , and ir^^^. The simultaneous fits of the K~tt + tt + and 
K + 7t~tt~ mass spectra with cuts optimized for D ± — > K^K^-k^ are shown in 
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Fig. §. The fits yielded 23465 ± 183 and 28014 ± 201 signal events for the D+ 
and D~, respectively. Fitting the K^-k + -k + and K + -k^ti^ mass spectra with 
cuts optimized for D ± -> (jm*, K*(892)K ± , and TT T n ± n ± yielded 31084 ±254, 
23765 ±182, and 20105 ±156 signal events for the D + , respectively. The signal 
yields for the D~ were 37057±281, 28384±199, and 23744±169, respectively. 
We investigated different techniques for extracting the number of D + and D~ 
signal events from the data of Fig. |5|, with special attention given to the excess 
of events seen just below the D mass. In the most extreme case, the ratio of 
D + to D~ increased by 1.9%, which changed the measured asymmetry by 
0.9%. This variation is sufficiently small that we neglected it in calculating 
the 90% confidence level limits. 

Using the measured yields of the various decay modes, we have calculated the 
CP- violation asymmetry parameters Acp directly from Eqs. ([3]) and (|J). The 
central values have statistical errors of the order of 3-5% and are all consistent 
with zero asymmetry. These results for Acp and the 90% confidence level limits 
are summarized in Table |l|. 

As discussed above, in order to extract Aqp with no bias, Eq. @ must be 
satisfied. Since we are normalizing to the CF decay D — > K t tc ± tt ± , most 
sources of systematic error cancel. Detector asymmetries may lead to a fake 
CP asymmetry. However, we expect such effects to be momentum dependent. 
We investigated any momentum dependence in our efficiencies by comparing 
the ratio of ratios of the yields of D + to D~ for all decays, as they appear in 
the efficiency relation given above in Eq. in bins of the D ± momentum. We 
also investigated momentum dependence in the asymmetry parameter itself 
by measuring Aqp for the different decay modes in bins of momentum. 
Within our statistical errors, we saw no momentum dependence in either test. 
Thus, we do not assign any systematic errors for these sources. 

Various additional sources of systematic uncertainty were investigated by vary- 
ing the cuts for the different decay modes, the and if* (892) mass sideband 
regions, and the fitting procedure, and by performing reflection subtractions 
instead of excluding the reflection regions from the fits. For all decay modes, 
the total systematic variations in Acp were found to be small compared to 
the statistical errors (less than half) and were thus neglected. 

In conclusion, we have measured Aqp for D ± — ► K T K ± ir ± , (pir^, K*(892)K ± , 
and 7r T 7r ± 7r ± . We see no evidence of anomalous CP violation in SCS charm de- 
cay. The statistical errors on our Aqp measurements for the — > K T K ± n ± , 
(p^, and K*(892)K ± are more than a factor of two smaller than the only 
previous measurements 0, and we have made the first measurement of Acp 
for the decay D ± — > 7r T 7r ± 7r ± . 
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Fig. 1. The fits of K~K + ir + and K + K~ir~ mass distributions, as described in the 
text. The Df and K^^ir^ reflection region above 1.93 GeV was excluded from 
both fits. The dashed lines represent the fits extended into the excluded regions. 
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Fig. 2. The fits of 07r + and 4>ir~ mass distributions, as described in the text. The 
Df region above 1.93 GeV was excluded from both fits. The dashed lines represent 
the fits extended into the excluded regions. 
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Fig. 3. The fits of K*(892)°K + and K*(892)°K~ mass distributions, as described in 
the text. The Df and K^^ir^ 1 reflection region above 1.93 GeV was excluded from 
both fits. The dashed lines represent the fits extended into the excluded regions. 
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Fig. 4. The fits of it~it + it + and 7r + 7r _ 7r~ mass distributions, as described in the 
text. The i^ =F 7r :t 7r ± reflection region below 1.83 GeV was excluded from both fits. 
The dashed lines represent the fits extended into the excluded regions. 
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Fig. 5. The fits of the K~ir + ir + and K + ir~7r~ mass distributions for i^ =F 7r ± 7r =F 
candidates that survive the D ± — > K T K ± ir ± selection criteria, as described in the 
text. 
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Table 1 

Summary of our measured CP-violation asymmetry parameters Acp and the 90% 
confidence level limits for the singly Cabibbo-suppressed decays D + — > K~K + tt + , 
07T+, K*(892)°K + , and tt-it+tt+. 



Decay Mode 


A C p 


90% CL Limits (%) 


KKir 


-0.014 ±0.029 


-6.2 < A C p < 3.4 


(j)TT 


-0.028 ± 0.036 


-8.7 < A C p < 3.1 


K*(892)K 


-0.010 ±0.050 


-9.2 < A C p < 7.2 


TTTTTT 


-0.017 ±0.042 


-8.6 < A C p < 5.2 
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